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The generalized Chaplygin gas provides an interesting candidate for the present accelerated 
expansion of the universe. We explore a geometrical explanation for the generalized Chaplygin gas 
within the context of brane world theories where matter fields are confined to the brane by means 
of the action of a confining potential. We obtain the modified Friedmann equations, deceleration 
parameter and age of the universe in this scenario and show that they are consistent with the 
present observational data. 
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CN ■ 1 Introduction 

> , 

In recent times the generalized Chaplygin gas (gCg) model has been proposed as an alternative to both 
the cosmological constant and quintessence in explaining the accelerated expansion of our universe. 
The Chaplygin gas model describes a transition from a universe filled with dust-like matter to an 
accelerated expanding stage. The generalized Chaplygin gas model, introduced in [1] and elaborated 
in [2], is described by a perfect fluid obeying an exotic equation of state 



where B is a positive constant and < f3 < 1. The standard Chaplygin gas corresponds to /3 = 1. 
An attractive feature of the model is that it can naturally explain both dark energy and dark matter 
[3]. The reason is that the Chaplygin gas behaves as dust-like matter at the early stages of the 
evolution of the universe and as a cosmological constant at late times. The Chaplygin gas appears as 
an effective fluid associated with (i-branes [4, 5] and can also be derived from the Born-Infeld action 
[6]. An interesting range of models have been found to be consistent with the SNe la data [7], CMB 
experiments [8] and other observational data [9]. The cosmological implications of the Chaplygin gas 
model have been intensively investigated in recent literature [10, 11]. 

Over the past few years, models with extra dimensions have been proposed in which the standard 
fields are confined to a four-dimensional (4D) world viewed as a hypersurface (the brane) embedded 
in a higher dimensional space-time (the bulk) through which only gravity can propagate. The most 
popular model in the context of brane world theory is that proposed by Randall and Sundrum (RS). 
The so-called RSI model [12] proposes a mechanism to solve the hierarchy problem with two branes, 
while in the RSII model [13] a single brane with positive tension is considered where 4D Newtonian 
gravity is recovered at low energies even if the extra dimension is not compact. This mechanism 
provides us with an alternative to the compactification of extra dimensions. 
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The cosmological evolution of such a brane universe has been extensively investigated and effects 
such as a quadratic density term in the Friedmann equations have been found [14]- [16]. This term 
arises from the imposition of the Israel junction conditions which is a relationship between the extrinsic 
curvature and energy-momentum tensor of the brane and results from the singular behavior in the 
energy-momentum tensor. There has been concerns expressed over applying such junction conditions 
in that they may not be unique. Indeed, other forms of junction conditions exist so that different 
conditions may lead to different physical results [17]. Furthermore, these conditions cannot be used 
when more than one non-compact extra dimension is involved. Against this background, an interesting 
higher-dimensional model was introduced in [18] where particles are trapped on a 4-dimensional 
hypersurface by the action of a confining potential V. In [19], the dynamics of test particles confined 
to a brane by the action of such a potential at the classical and quantum levels were studied and 
effects of small perturbations along the extra dimensions investigated. Within the classical limits, test 
particles remain stable under small perturbations and the effects of the extra dimensions are not felt by 
them, making them undetectable in this way. The quantum fluctuations of the brane cause the mass 
of a test particle to become quantized and, interestingly, the Yang-Mills fields appear as quantum 
effects. Also, in [20] a confining potential formalism was used to confine particles to an arbitrary 
manifold in a higher dimensional Euclidean space. In [21], a brane world model was studied in which 
the matter is confined to the brane through the action of such a potential, rendering the use of any 
junction condition unnecessary and predicting a geometrical explanation for the accelerated expansion 
of the universe. In a related work [22], a brane scenario was studied where the m-dimensional bulk 
is endowed with a Gauss-Bonnet (GB) term and the localization of matter on the brane is again 
realized by means of a confining potential. It was shown that in the presence of the GB term, the 
universe accelerates faster than brane models without the GB term. The behavior of an anisotropic 
brane world with Bianchi type I and V geometry was studied in [23] along the same lines. 

In this paper, we follow [21] and consider an m-dimensional bulk space without imposing the Z2 
symmetry. As mentioned above, to localize matter on the brane, assumed to be thin, a confining po- 
tential is used rather than a delta-function in the energy-momentum tensor. The resulting Friedmann 
equations on the brane are modified by an extra term that may be interpreted as the generalized 
Chaplygin gas, providing a possible phenomenological explanation for the accelerated expansion of 
the universe. It should be emphasized that the model presented in this work is different from that 
introduced in [24, 25] in that the latter provides no account for the confinement of matter to the 
brane. 



2 The model 

In this section we present a brief review of the model proposed in [19, 21]. Consider the background 
manifold V4 isometrically embedded in a pseudo-Riemannian manifold V m by the map y : V4 —> V m 
such that 

GA B yp B u =9^, Q AB y A ^ = Q, g AB N^AT b B = g ab = ±1, (2) 

where Qab (<7/xi/) 1S the metric of the bulk (brane) space V m (V 4) in arbitrary coordinates, {y A } ({2^}) 
is the basis of the bulk (brane) and M A are (m — 4) normal unit vectors orthogonal to the brane. 
Perturbation of V4 in a sufficiently small neighborhood of the brane along an arbitrary transverse 
direction £ is given by 

z A (x^,n = y A + (^y) A , (3) 

where C represents the Lie derivative and £ a (a = 1,2, ...,m — 4) is a small parameter along M A 
parameterizing the extra noncompact dimensions. By taking £ orthogonal to the brane, we ensure 
gauge independency [19] and have perturbations of the embedding along a single orthogonal extra 
direction J\f a , giving the local coordinates of the perturbed brane as 
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In a similar manner we see that since the vectors M A depend only on the local coordinates x M , they do 
not propagate along the extra dimensions. The above assumptions lead to the embedding equations 
of the perturbed geometry 

g^ = G AB Z A ^Z B ^ 9fia = G AB Z A ^ B a , G A BM A a N B b = g ab . (5) 

If we set Af A a = S A , the metric of the bulk space can be written in the following matrix form 

Gab = ( 9tIU + A A , (6) 

\ Kb g a b) 

where 

= 9»v - 2t; a K^ a + t a t h g^K ma K vfih , (7) 
is the metric of the perturbed brane, so that 

K^a = -GABy A nN B a ;v, (8) 



represents the extrinsic curvature of the original brane (second fundamental form). We use the 
notation A^ c = ^Ay^ where 

represents the twisting vector fields (the normal fundamental form). Any fixed £ a signifies a new 
perturbed geometry, enabling us to define an extrinsic curvature similar to the original one by 



Aa cd = GabM\M B c = \cd, (9) 



K^ a = -G AB Z A ^N B a . u = K^a - C b (K^ a K\ b + A^A C bu ) . (10) 
Note that definitions (6) and (10) require 

^ fj.ua ~ -o"^r- I 11 ) 



In geometric language, the presence of gauge fields A^a tilts the embedded family of sub-manifolds 
with respect to normal vectors J\f A . According to our construction, the original brane is orthogonal 
to normal vectors N A . However, equation (5) shows that this is not true for the deformed geometry. 
Let us change the embedding coordinates and set 

X A = Z A -g ab N A A b ,. (12) 

The coordinates X A describe a new family of embedded manifolds whose members are always or- 
thogonal to J\f A . In this coordinates the embedding equations of the perturbed brane is similar to 
the original one, described by equation (2) so that y A is replaced by X A . This new embedding of 
the local coordinates are suitable for obtaining induced Einstein field equations on the brane. The 
extrinsic curvature of the perturbed brane then becomes 

K, va = -GabX A »N b v = K, ua - t h K» ia K\ h = (13) 

which is the generalized York's relation and shows how the extrinsic curvature propagates as a result 
of the propagation of the metric in the direction of extra dimensions. The components of the Riemann 
tensor of the bulk written in the embedding vielbein {X A , Af A }, lead to the Gauss-Codazzi equations 
[26] 

Rap^s = 2g ab K a [ ia K S y 3b + UabcdX a X b X^X b ', (14) 
2if Q [ 7C;< 5] = 2g ab A^ ac K S ] ab + TZabcdX A M b X^X B , (15) 
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where IZabcd and R a p^s are the Riemann tensors for the bulk and the perturbed brane respectively. 
By contracting the Gauss equation (14) on a and 7 and defining 

Qiiv = ~9 ab (Kl a K lvb - K a K^ b ) + X - (K a p a K a ? a - K a K a ) 9fiu , (16) 

which is an independently conserved quantity, that is Q^ v v = [24]. Using the decomposition of the 
Riemann tensor into the Weyl curvature, the Ricci tensor and the scalar curvature, we obtain the 4D 
Einstein equations as 

G^v = GabX^X^ + Qfiu - £-nv + ( m - 2) ^^-AB^a^b 9\iv 

a b (to -4) 

(m — 2) ' F ' (to — l)(m — 2) 

where = 9 ab C ABCD-N'a is tne electric part of the Weyl tensor Cabcd- Now, let us 

write the Einstein field equations in the bulk space 

G% + K^Qab = a* Sab, (18) 

this equati 

consisting of two parts 



where a* = 1 _^ . In this equation is the cosmological constant of the bulk space with Sab 



Sab = Tab + \vGab, (19) 

where Tab is the energy-momentum tensor of the matter confined to the brane through the action 
of the confining potential V. We require V to satisfy three general conditions: it should have a 
deep minimum on the non-perturbed brane, depends only on extra coordinates £ a and the gauge 
group representing the subgroup of the isometry group of the bulk space should be preserved by it. 
Although the explicit form of such a potential is of no consequence in the present work, for the sake 
of completeness it would be useful to mention a result presented in [19] where a confining potential 
whose role, as the name suggests, is to trap test particles on the brane was obtained by invoking the 
above properties and assuming that the brane is located at £ a = 0, with the result 

V(6 = \w 2 g ab ei\ (20) 

where w is a constant much larger than the scale of curvature of the brane. This potential is clearly of 
the harmonic oscillator type, forcing a test particle leaving the brane along the extra dimension back 
to its initial position on the brane which coincides with the location of minimum of the potential. It is 
conceivable that other equivalent potentials having different forms but possessing the same properties 
and playing the same role can be envisaged. 
Using Einstein equations (18), we obtain 

n AB = —, °* x QabS + . 2 n M b) G A B + a* Sab, (21) 
(to — 2) (to — 2) 



and 



U= ( a *S -mA^). (22) 

TTXi — 



Substituting IZab and 1Z from the above into equation (17), we obtain 

r n 1 ( m ~ 3 ) ».*„abc n 1 2a * c (m - 4)(m - 3) n 

(to — 2) (to — 2) (to — 1)(to — 2) 



+ ^A<»V. (23) 
(to — 1) 
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As was mentioned in the introduction, localization of matter on the brane is realized in this model 
by the action of a confining potential. This can simply be done by taking 

2a* 

«V = ( m _ 2) Ta "" T ^ a = °' Tab = °' ^ 

where a is the scale of energy on the brane. Now, the induced Einstein field equations on the original 
brane can be written as 

G^u = ar^ v 2 ^ _ j-j OLTg^ v - Ag^ u + - £^ u , (25) 

where A = — (^Z^ A^ and Q^ u is an independently conserved quantity which may be considered as 
an energy- momentum tensor of the generalized Chaplygin gas (gCg). This matter is a candidate for 
a unified dark matter-energy scenario which is parameterized by an equation of state p = — , where 

B and are arbitrary constants. The speed of sound v 2 = — ^ in the gCg is defined as v 2 = — —. 
Ultimately, we have three different types of 'matter' on the right hand side of equation (25), namely, 
ordinary confined conserved matter represented by t^ u , the matter represented by Q^ u which will be 
discussed later and finally, the Weyl matter represented by £^ v . 

The geometrical approach considered here is based on three basic postulates, namely, the confine- 
ment of the standard gauge interactions to the brane, the existence of quantum gravity in the bulk 
and finally, the embedding of the brane world. All other model dependent properties such as warped 
metric, mirror symmetries, radion or extra scalar fields, fine tuning parameters like the tension of the 
brane and the choice of a junction condition are left out as much as possible in our calculations [25]. 
In the next section, we discuss the cosmological implications of our model. As was mentioned above, 
the results do not depend on the precise shape of the potential. 

3 Friedmann equations, deceleration parameter and age of the uni- 
verse 

In what follows we will analyze the influence of the extrinsic curvature terms on a FRW universe, 
regarded as a brane embedded in a 5-dimensional bulk with constant curvature (£^ u = 0). The FRW 
line element is given by 



ds 2 = -dt 2 + a(tf 



dr 2 9 
+ r 2 



(d6 2 + sin 2 9dip 2 ) 



(26) 



1 — kr 2 

The confined source is the perfect fluid given in co-moving coordinates by 

V = (P + P) u ^ + PflV«" u m = ~ S °^ P = (7 - l )P- ( 27 ) 
The constant curvature bulk is characterized by the Riemann tensor 

TZ-abcd = K(GacGbd - QadQbc)-, (28) 

where Gab denotes the bulk metric components in arbitrary coordinates and fc* is either zero for 
the flat bulk, or proportional to a positive or negative bulk cosmological constant respectively, corre- 
sponding to two possible signatures (4, 1) for the dS§ bulk and (3, 2) for the AdS§ bulk. We take, in 
the embedding equations, g 55 = e = ±1. With this assumption the Gauss-Codazzi equations reduce 
to 

Ra/3jS = -{KajK/35 - K a5 K (3y) + K(g a ~ t gf35 ~ 9aS9f3-y), (29) 
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K a{j3n] = 0. (30) 

Also the effective field equations derived in the previous section with constant curvature bulk can be 
written as 

= olt^ - Xg^ + Q^. (31) 

Here, A is the effective cosmological constant in four dimensions with Q^ v being a completely geo- 
metrical quantity given by 



£ 



{KK^ - K m K2) + X - (K ap K a P - K 2 ) g^ 



2 

where K = g^ v K^ u . Using the York's relation 



(32) 



K^a - 2 d ^ a , (33) 

we realize that in a diagonal metric, K^ va is diagonal. After separating the spatial components, the 
Codazzi equations reduce to (here a,/3,j,a = 1,2,3) 

K^ a + = K Vr + K i a r%, (34) 

K^ + = lK\ a . (35) 

The first equation gives K l la a = for a ^ 1, since K\ a does not depend on the spatial coordinates. 
Repeating the same procedure for a, 7 = i, % = 2, 3, we obtain K 2 2aa = for a ^ 2 and K 3 3a a = 
for <T^3. Assuming K\ a = K 2 2a = K 3 3a = b a (t), where b a (t) are arbitrary functions of t, the second 
equation gives 

Kooa = ~ (^f + b)j . (36) 

For /j,, v = 1, 2, 3 we obtain 

Kpva = Kg^. (37) 

Assuming that the functions b a are equal and denoting b a = b, h = | and H = f, we find from 
equation (32) that 



(38) 



3b 2 ( 2h \ 3b 
Qixv = — ~[Th + 1 ) 9 ^' M>^ = 1 > 2 ,3, <5oo = — • (39) 

As we noted before, is an independently conserved quantity, suggesting the same behavior as that 
of an energy-momentum tensor of an uncoupled non-conventional energy source. With this analogy 
we take the gCg model as a practical example and define Q^ v as an exotic fluid and write 

1 B 

Qfiu = - [(Pch + Pch)UnU u + Pch9nu] , Pch = p, (40) 
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where B is a positive constant and < [5 < 1 [2]. Comparing Q^v^p^v = 1,2,3 and Qoo from 
equation (40) with the components of and Qoo given by equation (39), we obtain 



Pch 



3^ / 2h_ 



Pch 



362 

ae 



(41) 



Equation (40) was chosen in accordance with the weak-energy condition corresponding to a positive 
energy density and negative pressure with e = +1. Use of the above equations leads to an equation 
for b(t) 



( 2ab N 
I 3d6 



B 



2* 

ae 



-(1+/3) 



(42) 



for which the solution is 



/ae \ 2 

m = ( T ) 



B + 



C 



a 3(l+/3) 



2(l+/3) 



(43) 



where C is an integration constant. Using equation (41) and this solution, the energy density of the 
gCg becomes 



Pch — Pch 



B s + 



a 3(l+/3) 



(44) 



where p c h is the gCg density at the present time, B s = Bp c j^ + ^ is a dimensionless quantity related 

to the speed of sound for the gCg today, v 2 s = (3Bp ch ^ + ^ and C = p^ h + ^ — B. From equation (44), 
we see that for B s = the gCg behaves like matter, p c h ~ a' 3 , whereas for B s = 1 it behaves as a 
cosmological constant, p c h ~ constant and when < B s < 1, the model predicts a behavior as that 
between a matter phase in the past and a negative dark energy regime at late times. This particular 
behavior of the gCg inspired some authors to propose a unified scheme for the cosmological "dark 
sector," an interesting idea which has been considered in many different contexts. 

The contracted Bianchi identities in the bulk space, G AB . A = 0, using equation (18), imply 



0. 



(45) 



Since the potential V has a minimum on the brane, the above conservation equation reduces to 



= 0. 



(46) 



For a perfect fluid as a confined matter source, the time evolution of the energy density of the matter 
is given by 



p a 



-37 



(47) 



Taking equation (31) and using the geometrical energy density for Q^u, the modified Friedmann 
equations on the brane become 



dV 



k A 



a 



~ ) + ~2 = q + oP0 a 37 + ^-Pch 
a J a z 3 3 3e 



a 



B s + 



(l-B s 



a 3(l+/3) 



(TT/3) 



(48) 



a A a 
a 3 o 



-37 



a 



(37 " 2) + -Pch 



B s + 



a: 



3(1+/?) 



-13 
(1+/3) 



B x 



(1-B S 



2a 3 (i+/3) 



(49) 
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Now, we rewrite the Friedmann equation (48) in terms of the redshift z = ^ — 1 and of the observa- 
tional parameter S7, as 



a 



h 2 = Hi n k (i + z f + n x + n m (i + zf^ + + (i - s s )(i + z) 3 ( 1+ « 



i 

i+0 



(50) 



where i?o is the present value of the Hubble parameter, fi m , and are respectively, the confined 
matter, cosmological constant and spatial curvature relative density parameters where the latter is 
associated with the geometrical dark energy through VL c h = ~^r- The deceleration parameter q 

defined by q = — p becomes 



n x + ^(l + z) 3 -^f Bs + (l-B s )(l + zf^ Bs -Ul-B s )(l + zf^ 



n k (i + z f + n x + n m (i + zf + ^f [b s + (i - b s )(i + *) w)] 



ir+sy 



(51) 



We consider £l\ = and show that, within the context of the present model, the extrinsic curvature 
can be used to account for the accelerated expansion of the universe. Figure 1 shows the behavior 
of the deceleration parameter as a function of the redshift for selected values of $7 m , B s and (3 for 
signature s = +1. Note that this behavior is much dependent on the range of the values that B s can 
take and is insensitive to the parameter (3. For having an accelerating universe with £l m = 0.3, the 
value of B s should lie in the range 0.5 < B s < 1. Also the value of B s determines the acceleration 
redshift z a . For B s = 1 and B s = 0.6 with Q m = 0.3, the accelerated expansion begins at z a = 0.67 
and z a = 0.07 respectively. 

Had we taken the bulk signature to be (3,2) or e = —1, equation (41) would have represented a 
fluid with negative energy and positive pressure, causing an unexpected behavior of the expansion of 
the universe. In order to better visualize this behavior, the deceleration parameter has been plotted 
as a function of the redshift for selected values of Q m and Q c h in figure 1. 

The age of the universe in FRW models is given by 



H Jo 



dx 



_, i > 

2 



(52) 



where Hq 1 = 9.8 x 10 9 h _1 years and the dimensionless parameter h, according to modern data, is 
about 0.7. Hence, in the flat matter dominated universe with fltotai = 1 the age of the universe would 
be only 9.3 Gyr, whereas the oldest globular clusters yield an age of about 12.5 with an uncertainly 
of 1.5 Gyr [27]. We find the age of the universe for our model by direct integration of the Friedmann 
equation (48) 

if 1 < 53 > 

Hq JO 



r 



,X-W + x 2Q ch 



R I VS.) 
-°s T x 3(l+/3) 



1 

1+/3 



1 ) 
2 



where 



Q,L. — — 



^0 



ZHf 



^ch — 



ml 



(54) 



For a flat, matter dominated universe with ft m = 0.3 and B s = 1 this leads to a prediction for the 
age of the universe of about 13.2 Gyr. It shows that the age of the universe in our model is longer 
than the FRW model. In figure 2 we have plotted the dimensionless age parameter Hoto as a function 
of £l m for some selected values of B s and f3. This behavior shows that for a fixed value of Q m the 
predicted age of the universe is longer for larger values of B s . Note also that the age parameter Hoto 
is an almost insensitive function of parameter (3 but that it depends strongly on variations of B s . 
This means that age considerations will be much more efficient for imposing constraints on the speed 
of sound, B s , than the values of the parameter (3 [28]. 
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q(z) q( z ) 




■/ 

-0.75,' 



-l' 2 4 6 8 10 Z 

Figure 1: Left, deceleration parameter as a function of the redshift for Q, m — 0.3 and B s = 1 (solid line), Q. m — 0.3 
and B s — 0.6 (dashed line), Q. m — and B s = 0.92 (dot-dashed line), Q, m — 0.1 and B s = 0.9 (dot-dot-dashed line) for 
signature e — +1 with (3=1 and right, the same parameter for signature e — — 1. 

Hoto 



2 

1 . 75 
1.5 
1 .25 
1 

. 75 
0.5 
0.25 




0.2 0.4 0.6 0.8 1 



Figure 2: Hoto as a function of fi m for B s = 1 (solid line), B s = 0.95 (dashed line), B s — 0.8 (dot-dashed line) and 
= 1. 

Finally, it would be appropriate to compare the predicted age of the universe in our model to 
the Randall-Sundrum brane model where the effects of the brane parameters and dark energy on the 
age of the universe have been studied. It has been shown that the effect of the quadratic term p 2 , 
resulting from the imposition of the Israel junction conditions on the energy density term is to lower 
significantly the age of the universe. This term effectively contributes as a positive pressure, making 
brane models less accelerating. This problem can be avoided if we accept the existence of dark energy 
with p = — |p (phantom matter) on the brane since it has a very strong influence on increasing the 
age of the universe [29] . 

4 Conclusions 

In this paper, we have studied a brane world model in which the matter is confined to the brane 
through the action of a confining potential, rendering the use of any junction condition redundant. 
We have shown that in a FRW universe embedded in a constant curvature dS$ bulk the accelerated 
expansion of the universe can be explained in a purely geometrical fashion based on the extrinsic 
curvature. We have established a correspondence with the phenomenological gCg dark energy model 
and have extended the predictions of geometrical matter in the more general case where the relation 
between p and p is not linear. Finally, we have found that the age of the universe in our model is 
longer than that predicted by the FRW models and hence more in line with the present observational 
data. 
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